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ABSTRACT 

We consider the statistics of pulsar binaries with white dwarf compan- 
ions (NS— WD). Using the statistical analysis method developed by Kim et al. 
(2003) we calculate the Galactic coalescence rate of NS— WD binaries due to 
gravitational-wave emission. We find that the most likely values for the to- 
tal Galactic coalescence rate [TZtot) of NS— WD binaries lie in the range 0.2-10 
Myr~^ depending on different assumed pulsar population models. For our ref- 
erence model, we obtain T^tot = 4.1ll|5g Myr~^ at a 68% statistical confidence 
level. These rate estimates are not corrected for pulsar beaming and as such 
they are found to be about a factor of 20 smaller than the Galactic coalescence 
rate estimates for double neutron star systems. Based on our rate estimates, we 
calculate the gravitational- wave background due to coalescing NS— WD binaries 
out to extragalactic distances within the frequency band of the Laser Interfer- 
ometer Space Antenna. We find the contribution from NS— WD binaries to the 
gravitational-wave background to be negligible. 

Subject headings: binaries: close-gravitational waves-stars: neutron-white dwarfs 

1. INTRODUCTION 



The observed properties of double neutron star (DNS) systems along with models of 
pulsar survey selection effects have been used for many years in order to estimate the coales- 
cence rate of DNS systems due to the emission of gravitational radiation (Narayan, Piran & 
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Shemi 1991; Phinney 1991; Curran & Lorimer 1995). In Kim, Kalogera, & Lorimer (2003; 
hereafter paper I), we presented a novel method to calculate the probability distribution of 
the coalescence rate estimates for pulsar binaries (see also Kalogera et al. (2004)). In paper I, 
we applied this method to Galactic DNS systems. Having a probability distribution at hand 
allowed us to calculate the most likely value for the Galactic DNS coalescence rate as well 
as statistical confidence limits associated with it. These were then used to calculate the ex- 
pected DNS inspiral detection rate for Laser Interferometer Gravitational- Wave Observatory 
(Abramovici et al. 1992). In this second paper, we extend our study to NS— WD binaries 
that are relevant to the future NASA/ESA mission Laser Interferometer Space Antenna 
(LISA; Bender et al. 1998) . 

There are now more than 40 neutron star-white dwarf (NS— WD) binary systems known 
in the Galactic disk (see e.g. Lorimer (2001) for a review). Here, we consider the subset 
of NS— WD binaries which will coalesce due to gravitational- wave (GW) emission within a 
Hubble time. There are currently three such coalescing binaries known: PSR J0751+1807 
(Lundgren, Zcpka, & Cordes 1995), PSR J1757-5322 (Edwards & Bailes 2001), and 
PSR Jl 141-6545 (Kaspi et al. 2000; Bailes et al. 2003). We calculate the Galactic co- 
alescence rate of NS— WD binaries based on their observed properties using the method in- 
troduced in paper I. The GW frequencies emitted by these systems fall within the ~ 0.1 — 100 
mHz frequency band of the LISA. Using our rate estimates we calculate the GW amphtude 
due to the NS— WD binaries out to cosmological distances and compare it to the sensitivity 
curve of LISA (Larson, Hiscock, & HeUings 2000) as well as the Galactic confusion noise 
estimates from white dwarf binaries (Bender & Hils 1990, 1997; Nelemans, Yungelson, & 
Portegies Zwart 2001; Schneider et al. 2001). 

The organization of the rest of this paper is as follows. In §2, we consider the hfetimes 
of NS— WD binaries and summarize the techniques we use to calculate their coalescence rate. 
The results of these calculations are presented in §3. In §4 we use our rate results to calculate 
the expected GW backgound produced by NS— WD binaries Finally, in §5, we discuss the 
implications of our results. 



2. COALESCING NS-WD BINARIES 

In general, the coalescence rate of a binary system containing an observable radio pulsar 
is defined by 

7^ = ^x^, (1) 

Tlife 
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where A^psr^ is the estimated number of pulsars in our Galaxy with pulse profiles and 
orbital characteristics similar to those of the known systems, /b is a correction factor for 
pulsar beaming and rnfc is the lifetime of the binary system. In the following subsections, 
we calculate the total lifetime of a pulsar binary and derive the probability density function 
(PDF) of the Galactic coalescence rate, -P(7^), for NS— WD binaries. 
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2.1. Lifetime of a NS-WD binary 

In Table 1, we summarize the observational properties and relevant lifetimes for the 3 
pulsar systems considered in this work. We define the lifetime of a coalescing pulsar binary 
Tiifc to be the sum of the current age of the observable pulsar and the remaining lifetime 
of the system. Assuming the pulsar spins down from an initial period Pq to the currently 
observed value P (both in s) due to a non-decaying magnetic dipole radiation torque (see 
e.g. Manchester & Taylor (1977)), its current "spin-down" age 

(2) 

where P (in ss~^) is the observed period derivative. For young pulsars like Jl 141— 6545, 
it is usually assumed that Pq <^ P so that Tgd reduces to the familiar characteristic age 
Tc = P/{2P) ~ 1.5 Myr. For the older recycled pulsars, however, Arzoumanian, Cordes, 
& Wasserman (1999) pointed out that this assumption is usually not appropriate, since 
the weaker magnetic fields of these objects mean that their present spin periods are only 
moderately larger than the periods produced during accretion. Adopting the spin-up line 
from Arzoumanian, Cordes, & Wasserman (1999), we may write 

/ p y/^ 

^0 = T? s . (3) 

Using the above two equations we calculate Tsa for the two recycled pulsars J0751-I-1807 and 
J1757-5322 to be 6.7 and 4.9 Gyr respectively (see Table 1). 

The remaining lifetime of a pulsar binary is defined by the shorter of the merger time 
of the binary due to the emission of GWs, Tmrg, or the time that the pulsar will reach the 
"death hne", (Ruderman &: Sutherland 1975). For young pulsars like J1141— 6545 which 
have relatively short radio lifetimes, T(j<Tinrg- Recycled pulsars, on the other hand, have far 



^This is equivalent to the so-called 'scale factor' (Narayan 1987). 
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smaller spin-down rates than young pulsars so that it is hkely that close binaries containing a 
recycled pulsar will coalesce before the pulsar reaches the death-line (Tmrg<Td)- For circular 
orbits, the results of Peters (1964) calculations for the merger time of a binary system of two 
point masses mi and m2 with orbital period can be written as simply: 

where the reduced mass fi = 17111712/ (mi -|- 1712). For the eccentric binary J1141-6545, we use 
the more detailed calculations of Peters (1964) to calculate Tmrg- Most of observed coalescing 
NS— WD binaries as well as the DNS systems have Tmrg^ 10^~^ yr. 

Our understanding of pulsar emission is rather poor and therefore it is not clear how 
to calculate an accurate time associated with the termination of pulsar emission and hence 
Td- Here we assume the spin-down torque is dominated by magnetic-dipole radiation with 
no evolution of the magnetic field. The surface magnetic field of a neutron star, Bg, can be 
estimated from the current spin period P (s) and spin-down rate P (ss~^): 

B, = 3.2 X lO^^Ppy/^ G . (5) 

Chen & Ruderman (1993) comprehensively discussed the evolution of a pulsar period based 
on different magnetic field structures. Their results are consistent with previous studies 
(Ruderman & Sutherland 1975; van den Heuvel 1987). We adopt their case C (eq. (9) in 
their paper) according to which the radio emission terminates when the "death-period" 

= ( 1.4x10^0 ) ^ 

is reached. Assuming that the surface magnetic field remains constant, we can integrate eq. 
(5) to calculate the time for the pulsar period to reach P^. We find that 

For PSR J 1141— 6545, we use eqs. (6) and (7) to find the remaining observable lifetime 
Td ~ 104 Myr. This is significantly less than Tmrg for this binary system (~ 600 Myr). 
Including the modest contribution from the characteristic age of Jl 141— 6545, we take the 
observable lifetime of the binary system to be riife=rc-|-rd~105 Myr. We note in passing 
that Edwards & Bailes (2001) estimated the remaining hfetime of PSR J1141— 6545 to be 
only ~10 Myr. Although no details of their calculation were presented in their paper, they 
probably assumed some decay of the magnetic field which led to their lower value Td and 
hence Tufe. 
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In the cases of the recycled pulsars J0751+1807 and J1757— 5322, which have lower 
magnetic field strengths and hence longer radio lifetimes, both binaries will coalesce before 
the pulsars stop radiating (i.e. rini.g<^ tj), so we calculate their lifetime using Tiifc=rsd+Tmrg- 
The estimated lifetimes arc 14.3 Gyr (J0751+1807) and 12.7 Gyr (J1757-5322), about two 
orders of magnitude longer than for the young NS— WD J 114 1—6545. 

2.2. Probability Density Function of the Galactic NS— WD Coalescence Rate 

Estimates 

The basic strategy we use to calculate PiJZtot) is described in detail in paper 1. In brief, 
using a detailed Monte Garlo simulation, for each observed NS— WD binary, we determine 
the fraction of the population that is actually detectable by careful modeling of all large-scale 
pulsar surveys. We include the selection effects that reduce the detectability of short-period 
binary systems when integration times are significant in comparison. Our pulsar population 
model takes into account the distribution of pulsars in the Galaxy and their luminosity 
function. Treating each pulsar seperately, our simulations effectively probe the specific pulsar 
sub-populations with pulse and orbital characteristics similar to those of PSR J0751-I-1807, 
J1757-5322, or J1141-6545. 

Prom the simulations we obtain N^hs pulsars detected by the surveys out of a Galactic 
population of A^tot in each model. We calculate A^obs repeatedly for a fixed A^tot- As shown 
in detail in paper I, the distribution of A^'obs follows a Poisson function, P(A^obs;<-^obs>)- 
We calculate the best-fit value of <A^obs>, which is the mean number of observed pulsars in 
a given sample, for a given A'tot- Since we consider each observed pulsar separately, we set 
-^obs = 1- For example, one PSR J0751+1807 and no other pulsar similar to this (in terms 
of spin and orbital properties of the pulsar) have been observed. The likelihood of detecting 
one pulsar similar to the observed one from the given pulsar population with A^tot samples 
is simply P(l;<A^ot,s>). We vary A'tot and calculate P(l;<A^obs>) to determine the most 
probable value of A^tot- Also, we found <A^obs> is directly proportional to A^tot- We calculate 
a, which is the slope of the function <A^obs>= ccA^tot for each observed system for a pulsar 
population model. 

Then, using Bayes' theorem, we calculate P(<A^obs>) from the likehhood P(l;<A^obs>) 
and eventually calculate P{Tl) using a change of variables. We repeat the whole procedure 
for all three observed coalescing NS— WD binaries, and combine the three individual PDFs 
to obtain a total PDF of Galactic coalescence rate of NS-WD binaries, P{Tltot)- 

In paper I, we showed that a normalized PDF of the coalescence rate for an individual 
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pulsar binary system can be written as follows: 

where Q is a coefficient determined by properties of the i*^ pulsar: 

^ /b 



(8) 



(9) 



Here, the beaming correction factor /b is the inverse of the fraction of An sr covered by 
the pulsar radiation beam during each rotation. In the case of the two DNS systems, 
PSRs B1913+16 and B1534+12, Kalogera et al. (2001) adopted fb ^ 6 based on pulse 
profile and polarization measurements of two pulsars. The lack of such observations for the 
current sample of NS— WD binaries means that it is difficult to estimate reliable values of 
fb. Therefore, in this paper, we do not correct for pulsar beaming (i.e. fh — 1). As a result, 
all our values should be considered as lower limits. 

In paper I, we calculated P{Tltot) considering two observed DNSs systems (labeled by 
the subscripts 1 and 2) . We defined the total rate = Tli + TZ2 and showed that 



C1C2 

C2 — Ci 



/ \C2 — Ci 



(10) 



where Ci < €2- In Appendix A, we show that this can be extended for the current case of 
interest where we have three binary systems such that 7^+ = TZi + 7^2 + "^3- This leads to 



(11) 



{C2-C,f{C^-C^f{C^-C2f 
(C3 - C2)3e-^i^+ [ -2(-2Ci + C2 + C3) + 7^+[-Cl(C2 + C3) + {Cl + C2C3)]] 

+ (C3 - Ci)3e-^^^+ [2(Ci - 2C2 + C3) + 7^+[C2(C3 + C,) - {Cl + C3C1)]] 
+ (C2 - Ci)^e-^^^+ [-2(Ci + C2 - 2C3) + 7^+[-C3(Cl + C2) + {Cl + C1C2)]] 

where the coefficients C\ (?' = 1, 2, 3) arc defined by eq. (9) and C\ < C2 < C3. This result was 
already used in our recent rate estimation for DNS systems to include the newly discovered 
pulsar, J0737-3039 (Burgay et al. 2003; Kalogera et al. 2004). As before, the confidence limit 
(CL)^ and the lower and upper limits (T^l and T^u) of the coalescence rate estimates are 



^Strictly speaking, the phrase "confidence hmit" used in this paper should be interpreted as a confidence 
interval such that a true value of the Galactic coalescence rate would exist in a given range of rates. However, 
we keep the terminology confidence limit for consistency with paper I and Kalogera et al. (2004). 
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defined in the same way we described in paper I, i.e. 



/ P(7^+) dn+ = CL , 



(12) 



and 



p(7eL) = P(7eu) . 



(13) 



3. 



RESULTS 



In Fig. 1, we show the resulting P{TZtot) for NS— WD binaries afong with the individual 
PDFs for each observed coalescing binaries. The figure shown here is obtained from our 
reference model (model 6 in paper 1). As we found for the DNS systems in paper 1, P(7^tot) 
is highly peaked and dominated by a single object. In this case, PSR Jl 141— 6545 dominates 
the results by virtue of its short observable lifetime (riifc ~105 Myr). This is in spite of the 
fact that the estimated total number of binaries similar to PSR J075 1+1807 (iVo75i ~ 2900) 
is the largest among the observed systems. 

Wc summarize our results for different pulsar population models in Table 2. The model 
parameters are identical to those described in paper I. We note however, following Kalogera 
et al. (2004), that our reference model is now model 6 (Lmin = 0.3 mJy kpc^) rather than 
model 1 (Lmin = 1-0 mJy kpc^). This choice reflects the recent discoveries of faint pulsar 
with 1400-MHz radio luminosities below than 1.0 mJy kpc^ (Camilo 2003). The peak values 
of P{TZtot) lie in the range between ~ 0.2 — 10 Myr~^ where the reference model shows a peak 
around 4 Myr~^. For the reference model, the uncertainties in the rates, (defined by lZ\]/TZi) 
are estimated to be ~ 6, 27 and 62 at 68%, 95%, and 99% CL, respectively. Comparing this 
to results from Kalogera et al. (2004), wc find that the uncertainties at different CL of the 
coalescence rate of NS— WD binaries arc typically larger by factor of ~1.4 than those of the 
DNS systems. This result is robust for all models we consider. 

The correlations between the peak value of the total Galactic coalescence rate T^pcak 
and the model parameters (e.g. the cut-off luminosity Lmin and the power index p) seem 
to be similar to those of DNSs we observed in paper I. As we found for the DNS systems, 
T^peak values strongly depend on a pulsar luminosity function rather than a spatial distribu- 
tion of pulsars in the Galaxy, in other words, T^peak values rapidly increase as the fraction of 
faint pulsars increases. 
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4. GRAVITATIONAL WAVE BACKGROUND DUE TO 

NS-WD BINARIES 

Close binaries consisting of compact objects (e.g. NS— WD binaries) are suggested as 
important GW sources in a frequency range below 1 niHz. In this range, due to the large 
number of sources, LISA would not be able to resolve each source within a given frequency 
band. Hence the Galactic binaries are expected to estabhsh a confusion noise level (or 
"background") dominated by WD-WD binaries (Bender &: Hils 1990, 1997; Nelemans, 
Yungelson, & Portegies Zwart 2001; Schneider et al. 2001). In this work, we consider the 
contribution from NS— WD binaries to the predicted confusion noise level. Using our results 
from the previous section, we calculate the amplitude of GW signals from NS— WD binaries 
in the nearby Universe and compare it with the LISA sensitivity curve^. In this work, we 
assume that the three observed systems represent the whole population of NS-WD binaries 
in our Galaxy. 

We calculate the characteristic strain amplitude of GWs (he) from NS— WD binaries 
using the results given by Phinney (2001) for circular binaries. In general, binaries with an 
eccentricity e emit GWs at frequencies / = ni', i.e. the n-th harmonic of the orbital frequency 
1/. In the case of circular binaries, n — 2 due to the orbital symmetry and the quadrupole 
nature of GWs. The eccentricity of PSRs J0751+1807 and J1757-5322 are e ~ 10"^ and 
~ 10~^, respectively. Hence, it is safe to consider them as circular binaries. PSR J1141— 6545 
has an appreciable eccentricity (e = 0.17), but for simplicity, we consider only the n — 2 
harmonic as if it were a circular binary^. 

For an observation of length Tobs with a GW detector, the contribution from background 
sources (NS-WD binaries in this work) depends on the number of sources within the fre- 
quency resolution, A/ = 1/Tobs- Following Schneider et al. (2001), we define an effective 

1 /2 

GW amplitude /irms(/) = ^c(/)(A///) ' , where hc{f) is a characteristic strain amplitude. 
Phinney (2001) showed a simple analytic formula to calculate hc{f) for a population of 
inspiraling circular-orbit binaries with a given number density in the nearby Universe. We 



■^We use the online sensitivity curve generator to calculate the sensitivity curve of LISA 
(http: / / www.srl.caltech.edu/ shane/sensitivity/MakeCurve.html) . 

*We calculate the power distribution in various harmonics for this eccentricity based on the result of Peters 
& Mathews (1963). We note that the GW amplitude calculated for Jl 141— 6545 in this work corresponds 
to ^70% of the total power of the gravitational radiation emitted from this binary. The remaining power is 
contributed from the higher harmonics. 
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use eq. (16) in his paper to calculate /ic(/)^, and find 

;^.(,)^1.7xlO-(^) (^) (^) (^) , (14) 

where is the "chirp mass" of a NS— WD binary defined by 

(Mns + Mwd)V5 ' ^'^^ 

and A'o is the comoving number density of NS— WD, i.e. the number of sources per Mpc^. 

We calculate the GW amphtude of NS— WD binaries in a frequency range /min < / < 
/max- Estimated GW frequencies of three NS— WD binaries based on their current sep- 
arations are all less than ~ 0.1 mHz. In our calculation, however, we set the minimum 
frequency /min to be 1 mHz taking into account the fact that the confusion noise level is 
mainly dominated by Galactic WD- WD binaries at lower frequency range / < ImHz (Nele- 
mans, Yungelson, & Portegies Zwart 2001). The maximum GW frequency /max is calculated 
by /max = 2/Pb,min, whcrc Pb,min is the minimum orbital period of the binary at the WD 
Roche-lobe overflow. 

Following Eggleton (1983), wc calculate the minimum possible separation of the binary 
using amin = Rwd/i^l, where tl is the effective Roche lobe radius. Wc estimate the radius of 
a white dwarf companion -Rwd adopting the results given by Tout, Aarseth, & Pols (1997). 
We show the estimated mass of white dwarf companions in Table 1. Gonverting Omin to /max 
based on the Kepler's 3rd law, we find that 



where Mch = 1.44M0 is the Ghandrasekhar limit. We note that /max = /max(^NS? ^wd)- 
Based on eq. (16), we define three frequency regions: (a) /min < / < /max, 0751, (b) 
/max, 0751 < / < /max, 1757, and (c) /max, 1757 < / < /max, 1141- In the rcgiou (a), for ex- 
ample, all three observed NS— WD systems contribute to the GW background. However, for 
frequencies larger than /max, 0751, PSR J0751-|-1807-like populations have already reached 
the Roche lobe overflow and we can not apply eq. (14) to these systems. Therefore, in a 
frequency region (b), we consider J1141— 6545-like and J1757— 5322-like populations. Simi- 
larly, for the highest frequency range (region (c)), we consider the contribution to the GW 
signals from PSR Jl 141— 6545-like population only. 



1 1/2 

^We assume a case that the last term in Phinney's eq. (16), (^^^^^-§^4 — becomes unity. The 
calculation is not significantly affected by different assumptions on cosmological models and comoving number 
density functions of coalescing binaries (Phinney (2001)). 
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in order to calculate /irms(f), we need the chirp mass and present-day comoving 
number density of NS— WD binaries A^o- Following Farmer & Phinney (2002), we define the 
"flux-weighted" averaged chirp mass: 



where jFg„ is the GW flux (^gw oc f^'^/^/A^^/-^ J^penk) and A/'peak is the peak value of P(A^tot) 
of the each sub-population of NS— WD binaries in our Galaxy. The subscript i repre- 
sents each pulsar sub-population. For example, in region (b), <Ai>= (A/i 757^^1757^^^^ -|- 
^flMlMn4l^^^^)/{^fl757Mn57^^^^ + Mui A<ii4i ^°^^) • Bccausc AApeak IS a constant and inde- 
pendent of the GW frequency, it follows that <M.> is independent of frequency. As a result, 
the evolution of orbital characteristics and hence the GW frequency of a binary are solely 
determined by the inspiral process. This is true regardless of the inital distribution of orbital 
characteristics. 

We now calculate the comoving number density of NS— WD binaries 



where N[z)dz is the number of NS— WD binaries per unit comoving volume between rcdshift 
z and z + dz. Noting that the number density of NS— WD binaries is proportional to the 
total number of systems, we may write 



where e is the star formation rate density per unit comoving volume (p) normalized to the 
Galactic star formation rate (r) i.e. e = p/r (in Mpc~^). A/psR is the most likely value of 
the total number of pulsars for each frequency range, (e.g. A/psR = A/0751 -|- A/1757 + A/1141 
in region (a)). We derived P(A^tot) for individual systems in paper I. In a similar fashion to 
the coalescence rate estimation described in Appendix A, the combined PDF of A^tot can be 
calculated from individual PDFs of the observed systems. Then A/psr can be obtained from 
the peak value of the PDF we calculate for each frequency range we discussed above. 

Using the results of Cole et al. (2001), we flnd 



<M>= 



E-^gw,i EAApeak,i A^i^"^/' ' 



(17) 




(18) 



No = e A/'psR , 



(19) 




(20) 



where {a,b,c,d) — (0.0166,0.1848,1.9474,2.6316) are parameters which take into account 
dust-extinction corrections (see Cole et al. (2001) for further details). Assuming a Hubble 
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constant Hq = 65 km Mpc~^, we calculate the Galactic star formation rate density 
p(0) ~ 0.01 Mq yr~^ Mpc~^. Following Cappellaro, Evans, & Turatto (1999), assuming 
the Salpeter initial mass function, we convert the Galactic supernova type SNii_,_ib/c rate^ 
to the star formation rate finding r ~ 0.7 MqYT~^. Numerically integrating eq. (20) out to 
^max = 5, which is considered to be the onset of the galaxy formation (Schneider et al. 2001), 
we find e ~ 0.6 Mpc~^. The number density of NS— WD binaries then can be calculated 
by eq. (19) for a given A/psR for each model. 

In Fig. 2, we plot the GW amplitude h„j^is against the simulated LISA sensitivity curve 
calculated for a signal-to-noisc ratio S/N=l. All dotted lines correspond to the GW am- 
plitude calculated from the full set of pulsar population models we consider and the solid 
fine is the result from our reference model. The range of the GW amphtude for all mod- 
els spans about an order of magnitude. We find that the GW background amphtude from 
NS— WD binaries is about 1-2 orders of magnitude smaller than the expected sensitivity 
curve of LISA at GW frequencies larger than ImHz and it is unlikely that this population 
will be detected with LISA. In the lower frequency region (below ~ ImHz), the GW back- 
ground amplitude from NS— WD binaries increases as / decreases. However, the contribution 
from NS— WD binaries to the GW background noise level would still be less than ~ 10% of 
the GW amplitude from WD- WD binaries (dashed line in 2). We note that, however, we 
have not considered any beaming corrections, so the NS— WD background curves should be 
viewed as lower limits. This possibility is discussed briefly in the next section. 



5. DISCUSSION 

We have used detailed Monte Carlo simulations to calculate the Galactic coalescence 
rate of NS— WD binaries. From the reference model, the most probable value of T^tot is 
estimated to be 4.11^2!5g Myr~^ at a 68% statistical confldence limit. We flnd that the 
coalescence rate of NS— WD binaries is about factor of 20 smaller than those of DNS for all 
pulsar population models we consider. As mentioned above, we did not take into account 
any beaming correction for NS— WD binaries. If we assume a beaming fraction of pulsars 
in NS— WD binaries similar to that of pulsars found in DNS, /b ~ 6, then the discrepancy 
beween T^peak (DNS) and T^pcak (NS— WD) is signicantly reduced. As a simple estimate, 
if we assume /b,ii4i ~ 5, but keeping /b = 1 for the other two binaries, the estimated 



^Cappellaro, Evans, & Turatto (1999) assumed a Hubble constant Hq = 75 km s~^ Mpc~^. Since we 
adopt Ho = 65 km s~^ Mpc~^, we have multiplied their results by a factor (65/75)^. We also note that we 
consider Sbc-Sd type galaxies only, which would be relevant for active star-forming regions. 
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Galactic coalescence rate increases to 18.06l^2!74 Myr"^ at a 68% confidence limit. Hence 
the ratio between DNS and NS— WD coalescence rate decreases to about 5. Because the 
contribution from PSRs J0751 + 1807 and J1757-5322 is an order of magnitude smaller than 
that of Jl 141-6545, moderate values of beaming fraction for those recycled pulsars do not 
change the result significantly. 

Based on the number of sources of NS— WD binaries in our Galaxy, we estimate the 
effective GW amplitude from the cosmic population of these systems. We find that the GW 
background from NS— WD binaries is too weak to be detected by LISA for the nominal 
beaming correction. Only by adopting an unreasonably large beaming correction factor, 
/b > 10, could these systems be detectable by LISA in the mHz range. These results are 
in good agreement with an independent study by Cooray (2004) based on statistics of low 
mass X-ray binaries. 

We finally note that combining the results from paper I and this work can give us strong 
constraints on the population synthesis models. The preferred models, which are consistent 
with both T^-NS-WD and TZbns, can then be used for the estimation of the coalescence rate 
of neutron star - black hole binaries, which have not yet been observed. 

We thank A. Cooray for noticing the frequency resolution correction factor in the calcu- 
lation of the GW background. This work is partially supported by NSF grant PHY-0121420 
and a Packard Fellowship in Science and Engineering to VK. DRL is a University Research 
Fellow funded by the Royal Society. He is also grateful for support from the Theoretical 
Astrophysics Visitors' Fund at Northwestern University. 
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In paper 1, we derived expressions for P(7^tot) for one and two coalescing binaries. In 
this paper, and in our revised DNS coalescence rate estimates (Kalogera et al. 2004), we 
extend this PDF to the case of three systems. Following paper 1, we define a coefficient for 
each observed NS— WD: 



where A < B < C. Recall that a is the slope of the function <A''obs>= aA^tot and is 
determined for each pulsar population model for each NS— WD system. By definition, the 
total Galactic coalescence rate is the sum of all three observed systems: 



Redefining 7^+ = T^a + "^b, where T^a = T^i + 7^2 and TZ^ = TZs, we transform TZ^ and 71^ to 
new variables TZ+ and TZ- = TZ^, — TZ\^ 



A. Combined rate PDF for three binary systems 




(Al) 



7^+ = 7^l + 7^2 + 7^3 . 



(A2) 



p(7^+,7^_) = p(7^a,7^b) 



dTZ- 



dTl+ 



(A3) 



This preprint was prepared with the A AS IM^rpC macros v5.2. 
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Since both 7?.+ and TZ- are positive, —TZ- < TZ+ < -\-TZ-. 

The PDF of the total rate 7?.+ is obtained after integrating P{TZ+, 1Z-) over TZ-: 

p{n+) = / p(7^+, 7^_) dn_ = ]- I p{n,, n^) dn^ , 



where 



(A4) 



(A5) 



Here, P{TZg, — TZi+TZ2) is given by eq. (10) and we can rewrite the formula with appropriate 
coefficients defined earher: 



e + e-^^^ 



.B-AJ I ^\ ) \B-A 

The individual PDF P{TZh) = PiT^s) (eq. (8)) is also rewritten as follows: 



(A6) 



(A7) 



Replacing 7^a = (1/2) (7e+ + TZ-) and 7^b = (1/2) (7e+ - 7^_) in eq. (A4), the normalized 
P{'Tl+) can be obtained by integration. After some algebra, we find: 



p(7^+) 



{B - A)^{C - Ay{C - B)^ 
(C - Bfe-^'^+ -2{-2A + S + C) + 7^+[-A(S + C) + (^1^ + SC)]] 

2{A -2B + C) + n+[B{C + A)- {B^ + CA)] 



(A8) 



3-BTZ4 



+ {C-AYe 



+ {B- Afe-^'^+ \-2{A + B-2C)+ n+[-C{A + B) + (C^ + AB)] 
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Table 1. Observational properties of NS-WD binaries. Prom left to right, the columns 
indicate the pulsar name, spin period P, spin-down rate P, orbital period P^, most probable 
mass of the WD companion m^d, orbital eccentricity e, characteristic age Tc, spin-down age 
Tsd, GW merger timescale Tmrg, time to reach the death line r^, most probable number of 
NS-WD systems of this type in the Galaxy and references to this system in the literature. 



PSRs 


P 


P 








e 






Tnirg 


Td 


A/pSR 


References'' 




(ms) 


(s s-i) 




(hr) 


(Mq) 




(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 






J0751+1807 


3.479 


8.08x10" 


21 


6.315 


0.18 


< 10"4 


6.8 


6.7 


7.6 




2900 


1,2 


J1757-5322 


8.870 


2.78x10" 


20 


10.88 


0.67 


lO"'' 


5.1 


4.9 


7.8 




1200 


3 


J1141-6545 


393.9 


4.29x10" 


15 


4.744 


0.986 


0.172 


1.5xl0"3 




0.6 


0.104 


400 


4,5 



''The assumed NS mass and inclination angle i are 2.2M0 and 78° for J0751-M807 (Nice et al. (2004)) and 1.35Mq and 
60° for J1757-5322. For PSR J1141-6545, we adopt values given in Bailes et al. (2003). 

''References: (1) Lundgren, Zepka, & Cordes (1995); (2) Nice, Splaver, & Stairs (2004) (3) Edwards, & Bailes (2001) ; (4) 
Kaspi et al. (2000) ; (5) Bailes et al. (2003). 
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Table 2. Estimates for the Galactic coalescence rate (JZtot) of NS-WD binaries at various 

confidence limits for all models considered. 



Model" Tetot (Myr-i) 

peak'' 68%'= 95%"= 



1 1.23 

2 1.00 

3 1.32 

4 1.53 

5 1.19 

6 4.11 

7 1.73 

8 0.83 

9 0.43 

10 1.42 

11 0.72 

12 0.55 

13 0.40 

14 0.23 

15 10.03 

16 3.68 

17 2.45 

18 1.55 

19 0.72 

20 5.48 

21 1.96 

22 1.55 

23 1.13 

24 1.14 

25 1.21 

26 1.33 

27 1.46 



+1.57 


+3.97 


-0.77 


-1.04 


+1.27 


+3.20 


-0.62 


-0.83 


+1.69 


+4.27 


—0.83 


— 1.12 


+1.99 


+5.02 


—0.97 


— 1.31 


+1.53 


+3.87 


—0.75 


— 1.01 


+5.25 


+13.23 


-2.56 


—3.47 


+2.21 


+5.57 


-1.08 


-1.46 


+1.07 


+2.71 


-0.52 


-0.71 


+0.55 


+1.39 


-0.27 


-0.36 


+1.88 


+4.76 


—0.92 


— 1.23 


+0.96 


+2.42 


—0.47 


—0.62 


+0.72 


+1.82 


—0.35 


—0.47 


+0.53 


+1.33 


—0.26 


—0.34 


+0.30 


+0.76 


—0.15 


—0.20 


+12.46 


+31.38 


—6.09 


—8.34 


+4.57 


+11.53 


-2.25 


-3.07 


+3.06 


+7.71 


-1.50 


-2.04 


+ 1.94 


+4.90 


-0.95 


-1.29 


+0.91 


+2.29 


-0.44 


-0.60 


+6.69 


+16.86 


-3.29 


-4.52 


+2.37 


+5.96 


-1.17 


-1.61 


+1.97 


+4.96 


-0.96 


-1.3 


+1.47 


+3.72 


-0.71 


-0.97 


+1.49 


+3.76 


-0.73 


-0.98 


+1.56 


+3.95 


-0.76 


-1.03 


+1.72 


+4.35 


-0.84 


-1.13 


+1.90 


+4.79 


-0.93 


-1.25 
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'^Modcl No. (see paper I for details 
of model parameters.) 

'^Peak value of the probability 
density function. 

'^Confidence limits. 
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Galactic coalescence rate (Myr~*) 



Fig. 1. — The PDFs of the Galactic coalescence rate estimation in both a logarithmic and 
a linear scale (inset) are shown for the reference model. The sohd line represents P(7^tot)- 
Other curves are P(7^) for PSRs J1757-5322 (dot-dash), J0751+1807 (short-dash), and 
J1141— 6545 (long-and-short dash)-like populations, respectively. Dotted hues correspond to 
68%, 95%, and 99% confidence limits for PiUtot)- 
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0.0001 0.001 0.01 0.1 
GW frequency (Hz) 

Fig. 2. — The effective GW amplitude /irms for coalescing NS— WD binaries overlapped with 
the LISA sensitivity curve. The curve is produced with the assumption of S/N=l for 1 yr of 
integration. Dotted lines are results from all models we consider except the reference model, 
which is shown as a solid line. We also show the expected confusion noise from Galactic 
WD-WD binaries for comparison (dashed hne). 



